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ABSTRACT 


Regenerators have been proposed to increase thermal cycle 
efficiency of gas turbine power plants which are being considered 
for installation in certain classes of United States Navy ships. 

These regenerators are to preheat intake air and are currently designed 
as compact heat exchangers. The recommended heat exchanger configur- 
ations are either counterflow or crossflow. This analysis considers 
only the counterflow configuration. 

Some of the structural problems associated with rectangular, com- 
pact heat exchanger design are discussed and certain structural design 
parameters are determined analytically. A digital computer program is 
presented for the analysis of the pressure effects in an idealized, 
two dimensional, rectangular heat exchanger. This program is written 
for a matrix ánalysis of an idealized heat exchanger with a cross 
section analogous to a multi-story building frame. 

This investigation was conducted by Robert L. Corbett and David 
W. Stubbs during the period from December 1965 to May 1966, at the 


United States Naval Postgraduate School, Monterey, California. 
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li. INTRODUCTION 

Several current projects in industry and one at the United States 
Naval Postgraduate School are actively concerned with the performance 
of gas turbine regenerators. The impetus for these investigations stems 
from a proposal to equip certain types of United States Navy ships with 
a power plant which will reduce, significantly, the delay time for 
transition from cruising speed to flank speed. 

This desired flexibility in speed is planned to be achieved either 
by augmenting or replacing the current steam turbine installations with 
diesel or gas turbine power plants. In general, the efficiency of the 
diesel is superior to that of the pure gas turbine installations; but 
the advantage in weight and space for gas turbine engines is significant. 
These advantages, coupled with increased efficiency gained through the 
use of regenerators, make the gas turbine appealing for shipboard 
installations. Therefore the design, operation, and maintenance of 
gas turbine regenerators are important in the proposed plan for gas 
turbine shipboard installation. 

Regenerator or heat exchanger requirements for different types of 
ships may differ in size, shape, and weight. A basic heat exchanger 
may be considered to consist of various combinations of small rectangular 
modules stacked or joined together in a manner that allows the complete 
heat exhcanger to meet the heat transfer, geometry, and weight require- 
ments for a particular installation. The rectangular shape is a log- 
ical choice for stacking and combining of modules and is simple to 
fabricate. 

Present design requirements specify maximum pressures up to ten 


atmospheres, temperatures ranging from 500 to 1000 ΘΕ, and a service 


life ranging from 20,000 to 60,000 hours. This latter requirement is 
severe, since structural integrity problems have occurred in similar 
heat exchangers after operating for much shorter periods of time. 

The Bureau of Ships has suggested that the following areas be 
investigated: 

1. Gas side fouling and air side salt precipitation 

2. Thermal stress behavior 

3. Thermal fatigue criteria 

4. Materials evaluation 

5. Header design and flow distribution 

6. Fabrication and joining techniques 

7. Advanced heat exchanger designs 

8. Engine evaluation of advanced heat exchangers 
This list indicates some of the areas which require investigation in 
order to define the requisite design criteria and to determine the order 
of reliability to be specified for a compact heat exchanger. 

Preliminary stress investigations of typical rectangular heat 
exchangers have been accomplished to date. Few complete analyses of 
the complex state of stress within the heat exchanger core or supporting 
frame work have been made and those reported are, in general, sufficient 
only for the justification of the preliminary design. These result 
in heavier and more costly heat exchanger design, and possibly in 
reduced reliability. The authors believe that, through a careful 
analytical study to determine certain structural design parameters, 
heat exchanger design may be determined more readily. 

The stresses in different parts of a heat exchanger are due 


primarily to the following: [2] 
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1. The difference in static and dynamic fluid pressures and 

atmospheric pressure. 

2. The thermal expansion or contraction of component parts 

under: 
a. Steady state conditions 
b. Transient or thermal shock conditions. 

3. Method of fabrication. 

The values of these stresses are difficult to determine not only 
because of the complex structure but also since the value of pressure 
and temperature is a function of location and time. 

The purpose of this investigation is to develop a basic method 
to study the effect of pressure on the design. Some of the problems 
relating to structural integrity are discussed and an attempt is made 
to isolate certain design parameters. The effect of some such para- 
meters is qualitatively studied through the use of the high speed 
digital computer. 

The investigation reported in this paper was conducted by the 
authors during the period from December 1965 to May 1966 at the United 
States Naval Postgraduate School, Monterey, California. Sincere 
appreciation is extended to the faculty and staff of the Aeronautical 
Engineering Department for their assistance. The authors are deeply 
indebted to Professor C. M. Smith for his invaluable encouragement and 
assistance throughout the project. The authors also wish to express 
their thanks to Professor C. H. Kahr for introducing them to and in- 


structing them in structural analysis by matrix methods. 
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2% GENERAL 

Heat exchangers of the cylindrical or tube type design have been 
analyzed and have proven to be structurally satisfactory. On the other 
hand, compact-rectangular heat exchangers have not been analyzed exten- 
sively for their structural action. 

A number of analytical approaches to the structural problem were 
explored in this investigation of pressure effects. Some of the methods 
proved to be not entirely satisfactory and sevetal of these are outlined 
in Appendix IT. 

The method of analysis adopted was limited to the core and was 
independent of the flange. The core is important from a heat transfer 
point of view, and is undoubtedly the most costly. The flange is not 
included in this investigation but its effect on structural design is 
sufficiently important to warrant further consideration. 

In the analysis of the core, cross sections were considered to be 
similar to multi-story building frames with the closure bars forming 
main structural columns, the fins acting as partition walls or columns, 
and the interior plates acting as floors. Exploratory design parameters 
were determined by considering a rectangular frame without interior 
tie-bars or constraints. The effect of these parameters on the lattice- 


like frame composing the core cross section was then investigated. 
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Su ANALYS IS 

For the purpose of this investigation, a າ. core module 
for a counterflow heat exchanger was postulated. In this module the 
geometry and material properties were assumed to be constant over the 
length; the pressure was assumed to be constant in all channels. The 
cross section was divided into a pattern of small rectangular passages 
by the interior plates and fins. All end discontinuities resulting from 
attachment of the core to flanges were neglected. This idealized core 


module is presented in Figure l. 





IDEALIZED COUNTERFLOW CORE MODULE 


This module represents only a small part of the complete regenerator 
heat exchanger. Different configurations would have a number of these 
modules arranged in a manner suitable to the operational environment 
of the ship. 

Close inspection of the idealized core module reveals that the 
fundamental element involved in the cross section is an open rectangular 
frame. Many of these basic frames are contained in a module cross section 
and two or more modules are stacked together to form a complete heat 
exchanger. 

Since the fundamental element is an open rectangular frame, an 
analysis was performed on a simple frame in order to isolate parameters 
which may affect or define structural considerations for design. [1] 

In this analysis the frame was assumed to be loaded with a uniform 
internal pressure and to have two axes of symmetry as shown in Figure 2. 
While these assumptions are not strictly applicable to the arrangement 
of elemental frames in the core module, they do facilitate the analysis 
and provide information as to which parameters affect the stresses in 


each member of the frame. 





FIGURE 2 
ELEMENTAL FRAME, TYPICAL 
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The two axes of symmetry which are shown in Figure 2 permit the 


analysis to be made for one quarter of the elemental frame; the required 


boundary conditions are defined through considerations of symmetry. 


The frame representing one quarter of the elemental frame is shown in 


Figure 3. | 1. 


Um jf 





FIGURE 3 
ONE QUARTER OF AN ELEMENTAL FRAME 


Considering member /\ of the frame in Figure 3, the definitive 


differential equation is; 


στ Тех (1) 
ақ) 


and the applicable boundary conditions are; 


a. Ji d Y, 


3 
ах? 


b. KE: = M, 
Хо 


С; qY, = O 
ах, ΧΟ 


= © 
χ.ο 





d. | = © 


х= % 
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dae 


Successive integration Qf Equation 1, with the application of 


boundary conditions yields: 


V X) - Qx, 
Мо = qx + М (2) 
Өк) E ax * Mya 

b d, Ji 


Now considering member A of the frame in Bigure 3, the definitive 


differential equation is: 


A. dh. = "9 (3) 
dx,* 


and the appropriate boundary conditions аге; 








д. Ja Pa = -H, = - OF 
dx? X270 > 
b. J2 а} = Ως” + M, 
dx X2=0 
c. AM] =Q%+MC 
ах іо 483 23, 
d. Ye = 0 
xX, 20 


Integration of Equation 3 and the introduction of the above 
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boundary conditions leads to the following expressions: 


Yep =a So. ae 


= 


Μω - 9% - аба + ຈ +”, (Ө 


So) 


да ов + @c + Әс? +М/ + с) 
6% t J> go 483, Jz 24 


From symmetry, the slope at the point X, = f/2 is equal to zero. This 
constraint allows the determination of the moment, Мі. 


M, = А |_2f?-3c# - (25) (5) 
om ++ c%G,) 


Assuming that the sides of the frame do not differ greatly in length 
or in bending stiffness, the maximum moment occurs at the corners of the 


frame. From Equation 2, the maximum moment can be expressed as: 


Mim = Qc J-3 4 2(&) - 3%) -ᾱ-5 (6) 
^u (Ye) + (GS) 


The maximum shear force on the horizontal side of the frame is: 


ка) мах = ຈ 


and the maximum shear force on the vertical side is: 


VOX) sx = га 


Equation 6 indicates that Ју? 


affect the choice of design from a structural point of view. Based on 


US f/c, and Q are parameters which 


the assumption that the pressure and side dimensions will be dictated 


by requirements for flow and heat transfer, f/c and Q may be fixed for 
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a given design for reasons other than structural. Thus the structural 
designer may have little control over the magnitude of the shear forces. 
The magnitude of the moments may be controlled by varying J, and ຽງ : 
Therefore the structural designer, by caréfully choosing the appropriate 
values of member thickness, may be able to optimize the design for the 
idealized frame element. 

In order to study the effect of varying member thicknesses, t1 


and t5, a two dimensional model of the heat exchanger module was in- 


vestigated. A diagram of this model is shown in Figure 4. This model 





FIGURE 4 


TWO DIMENSIONAL HEAT EXCHANGER MODEL 


shows that ty, t t_, and с. are analogous to the t 


а? ср 


Figure 3. The effect of varying ta» ty E and te within the model 


1 апа с, shown in 
chosen is not obvious and cannot be defined in terms of expressions 
derived for moment and shear for the simple frgme. 

In two dimensionalizing the model, the top and bottom plates, the 


interior plates, and the closure bars are in reality three dimensional 
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plates which have been idealized as two dimensional beams by assuming 
an elemental cross section of the core. The Plate БИРА! Rigidity, 
D, is used for each element to include this effect. 

Slope and deflection methods of analysis are adaptable to matrix 
methods which in turn are easily programmed for solution by the high 
speed digital computer. The computer program, as used, is included 
in Appendix I with a set of answers for the chosen model. This program 
Was prepared from considerations of the following. 

The stiffness method of matrix solution is a form of the slope 
and deflection method as proposed by M. L. Pei. [4] The basic element 
for this method consists of a beam under the action of end moments m 


1 


and m and end forces P, and P, as shown in Figure 5. 


CS wem 


FIGURE 5 
BASIC ELEMENT, SLOPE AND DEFLECTION ANALYSIS 
Let the beam undergo end rotations ©, and ©, as well as end displace- 
ments ô, and Θα. The positive directions of these displacements are 


as defined in Figure 6. 
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FIGURE 6 


BASIC ELEMENT DEFLECTIONS, SLOPE AND DEFLECTION ANALYSIS 
This beam then, has four degrees of freedom represented by the dis- 


placement matrix: 6 
| 


(ө) = (e 


5% 


Correspondingly, loads are represented by the load matrix: 


ч“ 


These loads and deflections are related by the stiffness influence co- 
efficient K which depends upon the cross section and the modulus of 
elasticity of the beam. These equations may be expressed conveniently 


in matrix notation as follows: 


{m} = [X] fe} 


For the single member example given above, this general equation 


becomes: 
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m, = КӨ + Кабг + Кз + Коф 
пг 2 Ко Ө, + K228: + Ke383 + Ko 94 
Ра: КӨ. + Kat Ot * Kss5s + Kay óy 
Py = Kai Dr + Kya 0, + Kusda + Кач 54 


Hence the "K" matrix is expressed as: 


Ки Кә Кз Кач 
[к] = Қи Қ Кгз Ка 


Қа! Κ.Σ Kas Καὶ 
kar Kya Куз К 


The loading (т) for the structure is known, so the rotations and 
deflections of the member ends, or more generally, the joints in the 


structure, may be determined from; 


е) (КТ [າ] 

The application of this method to a general two dimensional frame 
consists of computing the elemental "K' matrix for each element of the 
frame and then adding these elemental "K" matrices to form the total "К" 
matrix for the frame. Тһе "K" matrix is then inverted, as indicated by 
the above equation, to determine the frame joint deflections. 


Using the frame joint deflections the resulting element loads are 


[νὴ - [κ] (e) 


where [κ is the elemental matrix and (e) is the deflection of the 


joints in which the element terminates. These elemental loads are the 
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element internal loads. These internal loads describe the degree of 
relaxation from that of a fixed joint. The external loads or element 


reactions may then be determined from: 


Mpd Я M wreenas + M reaction 


These reactions are the desired results from which stresses may be cal- 
culated. 

The process of calculating loads in this way is very laborious, 
particularly when it is necessary to invert a large "КҚ" matrix. At the 
same time the process involved is purely a mechanical one which may be 


computerized. 
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4. RESULTS 
The computer was used to determine the effect of varying t ЕВ: 


апа 4 using the following variations: 


Basic Configuration Variations 

t, = .04 in .04 іп -- .25 іп 
ty = .04 in .04 in -- «25 іп 
tp = .04 in „02 in -- .06 іп 
tg = ¿Ol in Constant 

f = 0.75 in Constant 
ແ... Constant 


Maximum normal stresses and shear stresses for the variations shown 
were calculated. The shearing stress variations are functions of fin 
spacing and channel height, neither of which was varied since they were 
assumed to be dictated by thermodynamic considerations. The shearing 
stress variations were small when compared to the normal stress vari- 
ations. For this reason the variations in shearing stresses were assumed 
to be negligible for the parameters being investigated. Graphs showing 
the effect of the changes on shearing stresses are presented in Figures 
8 and 9. 

Maximum normal stress variations for varying closure bar thickness 
are plotted in Figure 10. Maximum stresses are computed from: 

Maximum normal stress - bending stress plus axial stress. 
As the thickness of the closure bar decreases; 

stress in the closure bar is increased rapidly, 

stress in the top plate is decreased rapidly, 

stresses in the interior plates and fins are decreased 


gradually. 
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Since the stresses in the closure bars increase rapidly, this method 
for reducing stresses in the model is considered marginal. 

Maximum normal stress variations for varying interior plate thick- 
ness are shown in Figure 11. This graph indicates that for an increas- 
ing interior plate thickness, the stresses in all members decrease 
rapidly. However, increasing the thickness of the interior plates 
rapidly reduces the heat transfer effectiveness between channels. Since 
heat transfer is the main purpose of the heat exchanger, this method 
of reducing the stresses must be used judiciously. 

Maximum normal stress variations for varying top plate thickness 
reduces the stresses in all members as long as the top plate is not too 
thick, see Figure 12. Beyond a certain top plate thickness the stresses 
in the top plate and the fins begin to increase. However, further thick- 
ness increase in the top plate is not necessarily detrimental as the closure 
bar stresses continue to fall off rapidly. At the point at which the 
top plate and the fin stresses reach a minimum and the closure bar stress 
reaches a local minimum, there is a significant shift in the location 
of maximum stress both in the closure bars and in the top plate. This 
maximum stress location shift results in a definite dip in the maximum 
stress values. Since variation of the thickness of the top plate has 
no basic undesirable effects, it is most advantageous to make initial 
stress reductions in the model by variation of top plate thickness. 

For any given heat exchanger this computerized method of solution 
is useful in the determination of the range of pressure stresses in 
the module. A computer program and detailed instructions for its use 


are presented in Appendix I. 
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5% CONCLUSIONS AND RECOMMENDATIONS 

The stresses due to internal pressure are a significant part of 
the stresses in a rectangular, compact heat exchanger. The value of 
these stresses for a pressure of 10 atmospheres may be on the order 
of sixty to ninety percent of design allowables. 

The value of stresses in the lattice-like frame which composes the 
heat exchanger core can be controlled through variations in the thick- 
ness of the individual elements involved. 

An optimum ratio of the sum of the stiffness of the top plate 
and interior plates to the stiffness of the closure bars is implied. 
This optimum should consider heat transfer as well as structural re- 
quirements. 

Consideration should be given to increasing the thickness of the 
top plate and/or the thickness of the interior plates in order to 
achieve the required stiffness for these members. This approach can 
lead to an optimum stress level and distribution for an internal pressure 
loading. 

The matrix analysis and the computer program used in this analysis 
should be expanded to provide solutions for three-dimensional problems 
with thermal as well as pressure loads. 

The compatibility of the flange and core is an important structural 
consideration. This should be thoroughly investigated in conjunction 


with further analysis of the core. 
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MAXIMUM SHEAR STRESS VERSUS TOP PLATE THICKNESS 


FIGURE 8 
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APPENDIX I 
COMPUTER PROGRAM '"FRAME'" DESCRIPTION AND USE 


Program Frame is a self-contained, Fortran 60, digital computer 





| program designed to provide joint deflections and/or rotation and member 
| end forces and/or moments for plane rigid frame type structures under 
general loading conditions. 

The geometry of the structure is restricted to a plane structure 
with all members parallel to either the X or Y axis in the usual 
cartesian coordinate system. 

A generalized flow diagram showing the basic mechanics of the 


program follows. 


Dimension the necessary matrices 
AK, À, C, THETA, AKS 


Read: NX - No. degrees kinematic freedom. 
NM - No. of frame members. 
EP - Error parameter, no. author is 
willing to accept as zero. 


Null A and C matrices. 


Read: M - No. representing moment deg. of freedom on left end of MBR. 
N - No. representing moment deg. of freedom on right end of MBR. 
K - No. representing shear deg. of freedom on left end of MBR. 
L - No. representing shear deg. of freedom on right end of MBR. 
EMOD - Elastic modulus of the member. 

ANER - Moment of inertia of the member. 


SPAN - Length of the member. 
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Compute elemental AK matrix. 
Store AK matrix in AKS matrix. 
Store M, N, K, L in AKS matrix. 
Build A matrix from AKM matrices. 


Transfer A matrix into C matrix. 


Read: NL - No. of loads on MBR. MBR not loaded 









Read: КЕМ - Concentrated or distributed load. 

W - Magnitude of external load. 
SPA 1 - Left distance from load to beam end. 
SPA 2 - Right distance from load to beam end. 








Member loaded 





1. 





Replace actual loads with fixed end member end loads and 
place in partitioned "C" matrix. 


2. Store fixed end member end loads in AKS matrix. 








Subroutine GAUSS 2. 


Partitioned "C" matrix inverted and multiplied yielding THETA 
matrix. 
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Take elemental "AK" matrix from AKS storage for one member, and 
multiply times DELTA (THETA for same member) to obtain internal 
moments and shears. 





Add internal moments and shears to those created by external 
loading (stored in AKS matrix) to obtain external end moments 
and shears. 












Print individual beam moments and shears. 


The following information is for those who wish to actually use 


the program in its present form. 

Properties such as the modules of elasticity (EMOD), the area 
moment of inertia (ANER), and the span may vary from one member to 
another; however, no provision is made for a single member with a 
variable moment of inertia (tapered member), except that an appro- 
priate average might be used if the operator so desires. 

Introduction of actual values of modulus of elasticity and moment 
of inertia will provide results, in appropriate units, of actual de- 
flections and/or rotations. Or, if desired, the relative magnitude of 
modulus of elasticity and moment of inertia, based on selected refer- 
ences may be used instead of using the actual values. In this case, 
the results will indicate the rotation or deflection times a factor 
which is equal to (Erof) х (1..р). In either case the correct values 
of member end forces and moments in appropriate units will result. 

The required input data is relatively straight forward and simple 
in that the first data card describes the overall structure and sub- 


sequent cards describe the individual members taken one at a time. 
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Care should be taken in compiling the data cards, since seemingly 
insignificant errors will produce entirely erroneous results. 

The program provides analysis for only one structure at a time. 
This structure may have only one compléte set of loads applied to it. 
The number of loads eomprising this one set, however, is essentially 
unlimited. Provision is made for application of concentrated loads 
and/or distributed loads to any member and in any combination. 

The program output format is self-explanatory except that THETA 
may be actual rotations or deflections if actual values of modulus of 
elasticity and moment of inertia were used, or they may be multiples 
of the reference values if such reference values were used. 

Following is the nomenclature of the program, a copy of the actual 
program in Fortran 60, a diagram for the model actually solved indicating 
member numbers and degrees of freedom (degrees of freedom indicated 
by arrows), data tabulation for the model as shown, a solution for 
deflections and rotations for the model, and the member end moments 
and shears for the model and data cards shown. 

Nomenclature - (Program Proper) 
AK - Stiffness matrix for an individual member. 
A - Matrix used for obtaining the stiffness matrix for the 
complete structure. 
C - Stiffness matrix for the complete structure. 
THETA - Matrix of deflections and/or rotations. 
AKS - Matrix used to store individual AK matrices and M, N, 
K, L, and end moments and forces due to external loads. 
DELTA - THETA matrix for a single member. 


F - Matrix of end forces and moments for a single member. 
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Required Input Data 


Data Card #1 Contains values of NX, NM, EP, punched in the 
format (2710, Е10.4). 

Data Card #2 Describes the first member to be considered. 

It contains values of M, N, K, L, EMOD, ANER, 
SPAN, punched in the format (412, 2E10.5, 
1F10.5). 

Data Card #3 Defines the number of loads (NL) applied to 
the member described in card #2. NL is punched 
in the format (12). If no load is applied to 
the member being considered, a Zero must be 
punched in column number 2. 

Data Card #4 Contains values of KFM, W, SPA 1, SPA 2, punched 
in the format (12 , 3F10.5). ΤΕ data card #3 
set NL equal to zero, this card (#4) must be 
omitted. The total number of #4 cards used at 
this point must equal NL from card #3. 

Subsequent Data Cards The input data card sequence is now 
repeated starting with card #2 above. Each 
member, in turn, must be completely described 
with its properties and loads before proceeding 
to the next member. 

Nomenclature - (Input Data) 

NX - Number of degrees of kinematic deficiency. 

NM - Number of members or spans comprising the structure. 

EP - Error parameter, minimum acceptable value of the deter- 

minant of the stiffness matrix. If determinant is less 


than E? the program will stop and print "matrix singular." 
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EMOD 


ANER 


SPAN 


KFM 


SPA 1 


SPA 2 


Basic Element 


The number assigned to the rotational freedom at the left 
end of a member. 

The number assigned to the rotational freedom at the right 
end of a member. 

The number assigned to the translational freedom at the 
left end of a member. 

The number assigned to the translational freedom at the 
right end of a member. 

Modulus of elasticity for a member. May also be used as 

a coefficient of some reference modulus. 

Area moment of inertia for a member. May also be used as 
a coefficient of some reference moment of inertia. 

The length or span for a member. 

The number of external loads (concentrated and/or uniform) 
applied to the single member being considered. 

Kind of load applied (1 = concentrated; 2 = uniform). 
Magnitude of the external load. 

Distance from left end of member to the location of a 
concentrated load. 

Distance from a concentrated load to the right end of a 


member. 


Note: Vertical members are rotated 90 degrees clockwise to determine the 


orientation indicated above. 
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The structure is sketched and all degrees of kinematic deficiency 
are numbered. 

The individual members are viewed one at a time, as shown on 
previous page, and M, N, K, and L for each member are set equal to the 
number which was used to designate the corresponding degree of freedom 
for that member. 


Load Parameters 


KFM = 4 
+W 
Γ SPA 4 T SPA? Ἴ 
| еч | 
КЕМ = 2 
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11 


20 


23 


26 


3615 
38 
40 


42 
43 


45 
46 


48 
49 


PROGRAM FRAME 
C RUN 110 


КЕАО 11» МХ»ММ» ЕР 


ҒОКМАТ 


(2110» 


Е10.4) 


DIMENSION AK(4s4),A(56»56),C(56,56), TH 


1DELTA(4)9F (4) 


DO 90 IO # 1, 


READ 20» M» N» 
2Е10.5» 


FORMAT 


(4129 


DO 23 121556 


DO 23 


J= 1956 


A(IsJ)= Oc 
00 26 Iz154 
DO 26 J=zl, 4 


AK(IsJ) = О 

АК (191) = (4 *XFEMODXANER)/SPAN 

АК (2+2) =АК(1»1) 

AK(1»2) = (2¢*EMOD*ANER) /SPAN 

AK( 193) * (64,*EMOD*XANER)/((SPAN**2) 
АК (293) = АК(1»3) 

AK(194) с -АК(1»3).: 

AK(294) = -AK(193) 

AK(393) = (12¢* EMOD*ANER) / ( SPAN*¥*3) 
AK(494) = AK(393) 

AK(394) = -AK(393) 


NM 


Ko Loe 


AK(291)=AK( 192) 
AK(391)#AK(193) 
AK(392)=AK(293) 
AK(491)=AK( 194) 
AK (492) =AK(294) 
AK (493)=AK( 394) 
AKS(10,17)=M 
AKS(10918)=N 
AKS(10»19)=K 
AKS(10,20)=L 


KKz1 


DO 3615 LL=194 
DO 3615 JJ=194 


AKS(109KK)=AK(LL» JJ) 


КК=КК+1 


IF (М) 
А(М»М) 


IF (NY 


AIMN) 
А(М»М) 
IF (K) 
A(M»K) 
А (КМ) 
IF (L) 
A(M»L) 
А((»М) 
IF (N 
А(М»М) 


103548538 
= AK(1»1) 
103, 42, 40 
= AK(1»2) 
= АК(1»2) 
103, 45, 43 
= AK(193) 
z AK(1»3) 
103, 48, 46 
= AK(1»4) 
= AK(1»4) 
103» 56 49 
= AK(292) 


ЕМОО» 
1F10.5) 


38° 


„4 


A A A e At 


ETA(56) sAKS(56924) » 


сет Фото тото сени nn nn юте — 


ΙΙ...“ a a =. rn ше 
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53 
54 


56 
DT 


59 
61 
63 


65 


681 
69 
πι 
73 

781 


76 
761 


78 
781 


80 
801 


82 
831 


85 
851 


87 
871 


89 
891 
90 


PRIM] 103,53, 51 

A(NsK) = AK(293) 

А(КзМ) = AK( 293) 

IF (L) 103, 56, 54 š 

A(NsL) = AK(2s4) 

AILIN) = AK(294) 

IF (K) 103, 61, 57 

A(KsK) = AK(393) 

IF (L) 103, 61» 59 

A(KsL) = AK(394) 

A(L»K) = AK(394) 

IF (L) 1039 639 62 

AlL»L) = AK(494) 

DO 65 Іг1%56 

DO 65 J=1956 

C(IsJ) * C(I»J) + Alls J) 

NPM=NX +1 

READ 681» NL 

FORMAT (12) 

IF (NL) 103, 90, 69 

00 90 KO = 1»NL 

READ 71» KFMs Ws SPAl» SPA2 

FORMAT (1I2» 3F10.5) 

IF (KFM-2) 73, 82, 103 

IF (M)103, 741» 74 E 

С(МӘМРМ) * (W*SPA1*SPA2**2)/(SPAN**2) + C(MsNPM) 
AKS(IO»s21)2-(W*SPA1*SPA2**2)/(SPAN**2) +AKS(10921) 
IF (N) 1039761976 | 
С(М,МРМ) * -(W*SPA2*SPA1**2)/(SPAN**2) * C(N5»NPM) 
AKS(IO$22)(W*SPA2*SPA1**2)/(SPAN**2) ғАҚ5(10,22) 
IF (K) 103,781,78 

С(КӘМРМ) - ((МҰ5РА2Ұ%Ұ2)%( 3,Ж5рРДІ + SPA2))/(SPAN**3) + C(K»NPM) 
АҚ5( І10»23)--((МЖ5РА2%%2)%( 3.%5рРА1%5РА2))/(5РАМЖҰҰ3) +AKS(109 23) 
ІҒ (() 103»801»80 

C(L3NPM) s ((W*SPA1**2)*(3,.*SPA2 +: SPA1))/(SPAN**3) + C(L»NPM) 
ο ο ο στι κ a a ee eet) eas +AKS( 10924) 
GO TO 90 

IF (M) 103582315823 

C{M»NPM) = (W*SPAN**2)/(12.) % С(МӘМРМ) 

AKS (10921) =—(W¥SPAN*®#2)/(12¢) +AKS(10»2]) 

ТЕ (N) 10358515985 

C(N»NPM) = -(W#SPAN**2)/(12.) + С(М,МРМ) 

АК5 (10, 22) = (МҰ5РАМҰҰ2)/(12.) +AKS(10922) 

IF (K)1039871987 

C(K»NPM) = (W*SPAN)/(20) + C(K NPM) 
AKS(10923)=—(W*SPAN) /(2e) +AKS( 10923) 

IF(L) 10398915389 

C(LÍNPM) = (WXSPAN)/(20):+ С( 9 МРМ) 
AKS(10924)=-(W*SPAN)/(20) +AKS(109,24) 

CONTINUE 

NP = 1 

CALE GAYUSS2 (NXy NPy EPs ECE»; THETA; Kł) 

GO TO (93, 99) » Kl 


93 
94 


321 


96 
97 


9072 
9074 
73230 


731 


976 


979 


9711 
9712 


9714 


9716 
9717 


9719 


9721 
9722 


9724 
9726 
9727 


9730 


9735 
9736 


99 


PRINT 94 | 

FORMAT(////*4H1NOTE le THETA IS JOINT ROTATION/DEFLECTION/) 
PRINT 321 

FORMAT (29H NOTE 2e ANSWER = EI X THETA////) 


DO 96 K=1»NX 

PRINT 979 Ks THETA(K) 

FORMAT (11HOEI X THETA 13»3H = F20.15) 

PRINT 9072 

FORMAT (////43HIMEMBER END MOMENTS AND SHEARS APPEAR BELOW/) 
PRINT 9074 

FORMAT(47H CLOCKWISE MOMENTS AND DOWN SHEARS ARE POSITIVE////) 
PRINT 730 : ' y н 
FORMAT(////43H 

PRINT 731 

FORMAT (43H MEMBER 
DO 9735 IM=],»NM 

KJ= 1 

DO 976 LL=1,4 

DO 976 JJ=1»4 
AK(LL+JJ)=AKS(IM9KJ) 
KJ=KJ +] 
M=AKS(IM917) 

IF (M) 103997119979 
DELTA(1)=THETA(M) " 
GO TO 9712 

DELTA(1)=08 

N=AKS(IM» 18) 

IF(N) 1039971699714 

DELTA(2)= THETA(N) 

GO TO 9717 

DELTA(2)=0. 

K=AKS(IM» 19) 

IF (К) 1039972199719 

DELTA(3) =THETA(K) 

GO TO 9722 $ 
DELTA(3)=08 N 
(-АК5(ІМ»20) 

IF (L)103,9726,9724 
DELTA(4)=THETA(L) 

GO TO 9727 

DELTA(4)=08 

DO 9730 KB=1 4 

Е(КВ) =0, 

DO 9730 JB=1»4 
F(KB)=F(KB) +AK(KBoJB) 
F(1)=F(1) +AKS(IM>21) 
F(2)=F(2) + AKS(IMs22) 
F(3)=F(3) + AKS(IMs23) 
F(4)=F(4) + AKS(IM»,24) 
PRINT 97369 IMoF(1)9F(2)9F (3) oF (4) 
FORMAT (149 F9.393F70.3) 

GO TO 105 

PRINT 100 


LEFT RIGHT LEFT RIGHT) 


MOMENT MOMENT SHEAR SHEAR//) 


* DELTA( JB) š 


&0 


Se a EE ——— 


Ч Tu. Ge ш ЖЕ —s sa 


100 FORMAT (16H MATRIX SINGULAR) 
GO TO 105 
103 PRINT 104 ; 
104 FORMAT (12H INPUT ERROR) : 
105 CONTINUE 
STOP 
END 
SUBROUTINE . GAUSS2 (NM. EP, A, X »KER) 
DIMENSION A(56956) 9X(56) 
NPM=N+M 
10 DO 34 L=1»N 
KP=0 
220.0 
| 00 12 KsL,N 
x ІҒ(2-АВ5Ғ(А(К»())211912912 
11 7-АВ5Ғ(А(К»0)) 
g KP=K 
| : 12 CONTINUE 
| IF(L-KP)13520,20 
| 13 DO 14 JsL,NPM 
Z=A(LoJ) 
A(LyJ)=A(KP sJ) 
14 А(КР».Ј) =2 
20 IF(ABSFIA(L+L))-EP)50>50930 
30 IF(L-N)315405»40 
21 LPI=L+] 
DO 34 K=LP1s»N 
IF(A(KsL)) 32934932 
32 RATIO=Al(KsL)/A(LoL) 
DO 33 J=LP1>NPM 
33 A(KyJ)=zA(KyJ)-RATIO*A(L +J) 
34 CONTINUE 
40 DO 43 Is1»N 
II=N+1-1 
DO 43 J=1»M 
JPN=J+N 
S=0.0 
IF(II-N)41543543 
4) IIPlsII*1 
DO 42 K=IIPl»N 
42  S=S+A(II>K)*X(K) 
43 X(II)= (А(11»УРМ)-5)/А(11911) 
KER=1 
RETURN 
50 KER = 2 
END 
END 
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APPENDIX IT 
OTHER METHODS OF ANALYSIS 

Prior to the decision to use the two-dimensional frame analysis 
for the internal pressure stresses a number of other methods were in- 
vestigated. Generally, the approach was one of searching for the correct 
design philosophy and attempting to adapt simple solutions to the prob- 
lem at hand. 

The methods outlined herein seemed to apply to the problem but in 
each case the results became too complex to provide useful design in- 
formation or the results required specific knowledge of some particular 
parameter which was unavailable. 

TOP PLATE ON AN ELASTIC FOUNDATION 

The first approach was to consider the top plate of the heat ex- 
changer to be fixed on the two sides adjoining the side closure bars, 
to be resting on an elastic foundation, with modulus K lb/in, along its 
span, and to be loaded with a uniform pressure load, Q lb/in. The elastic 
foundation was assumed to be the result of the arrangement of interior 
plates and fins. Assuming no change in geometry or other parameters 
with length, the plate was analyzed as a beam considering the stiffness 
to be EI/(1- v2) instead of simply EI. 


The expression for the moment at the end of this beam is: 


МЕ de et | 


2A? LSWRAD « SIN Ab 


where: 4 K(-v2) 
“: / ມຸ 
This result shows that the modulus of the elastic foundation, K, is of 


primary importance. [In truth, this modulus is a function of such para- 


meters as fin spacing, width of the heat exchanger, and thickness or 
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stiffness of the interior plates. No adequate analytical means of 
estimating the magnitude of the foundation modulus was found. It is 
recommended that K be determined by testing actual modules. 
DISCONTINUITIES BETWEEN HEADER AND CORE 

A number of attempts were made to determine the value of discon- 
tinuity stresses at the joint between the header or flange and core. 
In general, this method attempted to extend the classical solutions for 
cylindrical pressure vessels to include rectangular pressure vessels. 
This approach might be possible if the corners of the rectangular heat 
exchanger could be rounded off slightly such that a continuous mathematical 
function could be used to define the surface of the heat exchanger. The 
best approach with the square corners appears to be with the use of an 
elemental lengthwise strip taken at the location of the greatest dis- 


continuity in deformation between the core and the flange as shown below. 





If Q, and My could be determined, the elemental strip of the top 
plate of the core could be analyzed as a beam on an elastic foundation 
with end loads Q, and M,- Difficulties, in connection with finding Qo 
and Mọ» the magnitude of the elastic foundation modulus, whether the 
elemental strip of the core could be considered semi-infinite in length, 
and the correct method for determining the d between the header and 


core, were encountered. 
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FRAME ANALYSIS 

A number of lattice-type frames were analyzed and the analysis for 
the simple frame with two axis of symmetry did, in fact, provide the 
impetus for continuing with the complex frame analysis as used. This 
analysis provided fundamental information as to the effect of the various 
parameters involved in the core. 

Frames with various numbers of non-yielding tie bars, which repre- 
sented the interior plates, were analyzed. The solutions neglected the 
effect of fins. This approach was abandoned because of the inability 
to adequately account for the effect of the fins. 

SANDWICH CONSTRUCTION 

Considerable effort was expended in researching publications dealing 
with sandwich construction. Much is available on the subject in the 
literature but no reference was found which dealt with sandwich or honey- 


comb structures loaded with an internal pressure. 
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